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Abstract 

Two-dimensional spin-uncompensated momentum density distributions, Pg D (p)s, were recon- 
structed in magnetite at 12K and 300K from several measured directional magnetic Compton 
profiles. Mechanical de- twinning was used to overcome severe twinning in the single crystal sample 
below the Verwey transition. The reconstructed pg D (p) in the first Brillouin zone changes from 
being negative at 300 K to positive at 12 K. This result provides the first clear evidence that elec- 
trons with low momenta in the minority spin bands in magnetite are localized below the Verwey 
transition temperature. 
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I. INTRODUCTION 



Magnetite, Fe3C>4, is the oldest magnetic material recognized 2500 years ago and is one of 
the most fascinating materials in present-day solid state physics. It is because Fe 3 4 exhibits 
many interesting physical properties such as mixed valence and a metal-insulator (MI) 
transition known as the Verwey transition j^]. Furthermore, Fe3C>4 has potential industrial 
applications in magnetic multilayer devices jj| since it has full spin polarization with a 
high magnetic ordering temperature and its electronic structure at ambient conditions is 
predicted by band calculations to be half- metallic Q]. 

Above the MI transition temperature, T v , Fe3C>4 has an inverse spinel structure with 
cubic symmetry where eight Fe 3+ ions occupy the A sites with local tetragonal symmetry, 
while the B sites, which locally exhibit octahedral symmetry, are occupied by eight Fe 3+ 
and Fe 2+ ions in a simple ionic chemical formula. The results of resonant x-ray scattering 
{(J show that the Fe ions on the B sites are electronically equivalent for time scales lower 



Dy a structural change 
3]. Early explanations 



than 10~ 16 s. The MI transition occurs at T v ~ 120 K accompanied 
where the conductivity decreases sharply by two orders of magnitude 
of the MI transition invoked charge ordering on the Fe B sites below T v . However, NMR 
relaxation results indicate that the states of Fe ions on the B sites are mixed strongly even 
below T v [7jl. Moreover, recent resonant x-ray and inelastic neutron scattering studies show a 
fractional charge orde ring on the Fe B sites below T v |8|, [9|, 
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band calculations 
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12| , which is predicted by 



14j . Thus, some studies suggest that only the Fe B sites are involved 



in the change in valence across the Verwey transition [l5l. 1 1 61 ] while other investigations 



claim that the Fe A sites play also an important role |17l . |18| . Therefore, despite continuing 
experimental and theoretical efforts in the last 70 years, many questions regarding the MI 
transition still remain open. 

Magnetic Compton scattering using synchrotron radiation is an established technique 
for probing the spin- uncompensated momentum density distribution in a material with 
spontaneous magnetization [19|, |20J]. The total electron momentum density, p(p), is 

p(p) = pt(p) + Pi(p)> 

where p^(p) and pj,(p) give the momentum densities of electrons in the majority and minority 
spin bands, respectively, and p is the electron momentum. Owing to the presence of spin- 
dependent terms in the scattering cross- sect ion, the spin-uncompensated momentum density 



distribution, p s (p) = Pt(p) — P|(p)> can be measured using circularly polarized x-rays. It 
is obtained in the form of a projection of p s (p) on the x-ray scattering vector, the so-called 
magnetic Compton profile, defined as 



Here the z-axis is taken to be parallel to the x-ray scattering vector and p z is the projection 
of p on the z- axis. The integral of J m ag(Pz) over p z yields the spin- magnetic moment of the 
material. For electronic structure studies in a material with spontaneous magnetization, a 
measurement of p s (p) is very useful because p s (p) reflects the spin polarization of occupied 
bands in momentum space. It is possible to image p s (p) from a finite number of J ma ,g{Pz)s 
using a reconstruction technique. 

In this work we report magnetic Compton scattering studies of a single crystal of mag- 
netite. The amplitude of the two-dimensional spin-uncompensated momentum density dis- 
tribution is shown to develop negative excursions above the Verwey temperature. We con- 
clude on this basis that a substantial amount of t 2g minority electrons become itinerant on 
Fe B sites above T v . 

II. EXPERIMENTAL 

A high-quality single crystal of Fe3C>4 was synthesized using the floating zone melting 
method. The crystal was annealed at around 1200 °C in the controlled atmosphere of a 
mixture of CO and CO2 in order to homogenize the composition of the specimen and to 
release internal strain. We used a disk-like single crystal sample of approximately 10 mm 
diameter and 3 mm thickness with its axis parallel to the cubic [101] axis. 

Magnetic Compton scattering experiments were carried out at 12 K and 300 K on beam- 
line BL08W at SPring-8, Japan. Elliptically polarized x-rays emitted from an elliptical 
multipole wiggler were monochromatized to 174 keV by a bent Si(620) crystal. Energy 
spectra of Compton scattered x-rays from the sample with a scattering angle of 178.0° were 
measured using a 10-segmented Ge solid-state detector with external magnetic field of 20 
kOe. The estimated momentum resolution is 0.50 ± 0.01 a.u. full width at half-maximum. 

The directional J ma ,g(p z ) is extracted from the difference between two spectra taken under 
the same experimental conditions with alternating directions of magnetization of the sample, 
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aligned by an external magnetic field. The observed spectra were corrected for the energy- 
dependent scattering cross section, efficiency of the detector and absorption of the sample. 



III. RESULTS 

Since severe twinning occurs at T v in the single crystalline sample due to the structural 
change, de-twinning of the sample below T v was accomplished by simultaneous application 
of external magnetic field cooling and squeezing to extract directional magnetic profiles 
JmngiPz) at 12 K. An external magnetic field of 20 kOe was applied along the cubic [001] axis 
to establish the c-axis below T v . The use of a special Cu sample holder enabled us to squeeze 
the sample along the cubic [111] direction as a result of the larger Cu thermal contraction. 
Changes in the extracted J m!ig (p z ) at 12 K were observed below \p z \ ~ 2 a.u. resulting 
from this mechanical de-twinning. Thus the mechanical de-twinning is effective although 
we have not estimated the twin ratio in the sample at 12 K with and without mechanical 
de-twinning. Figure 1 shows the extracted J ma , g (p z ) along the three cubic principal axes at 
12 and 300 K. These profiles show a substantial directional dependence of the magnetic spin 
momentum density in good accord with the corresponding theoretical predictions discussed 
below. 

Figure 2 shows the extracted J ma , g (pz) at 12 K and 300 K between the cubic [001] and 
[101] axes. Note that the dip in J ma g(p z ) below p z ~ 1 a.u. observed from 9° to 27° at 
300 K disappears at 12 K. This result reveals that p s (p) changes across T v even though the 



saturation magnetic moment changes by less than 0.1 % at T v 21] . We have reconstructed 
the two-dimensional spin-uncompensated momentum density distribution, p^ D (p). For this 
purpose, we adopted the so-called direct Fourier-transform method to reconstruct p^ D (p) 



from several directional J m ag(p 2 )s [22j, [23]. In the reconstruction procedure, we assumed the 



cubic [010] axis to possess a four-fold rotational symmetry, even though Fes04 at 12 K has 



a lower than cubic symmetry 2J, |25j. This assumption was made as we measured J mag (p 2 )s 
only from the cubic [001] to the [101] axis due to the mechanical de-twinning using the 
special Cu sample holder. 

The reconstructed p^ D (p)s are shown in Figs. 3(a) and (b), which correspond to projec- 
tions of p s (p) on the cubic (010) plane. The reconstructed Ps D (P(oio)) looks like a volcano at 
12 and 300 K, where P(oio) represents the projection of p on the cubic (010) plane. At 300 
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K, the reconstructed /^(PCoio)) i n the crater region (below |p[ooi]| and |p[ioo]| ~ 0-4 a.u.) is 
negative, while highly positive values are found in the region of |p(oio) | ~ 1-7 a.u. with peaks 
at |P[ooi]| — |P[ioo]| ~ 1-2 a.u., where p[ooi] and P[ioo] represent the components of p along 
the cubic [001] and [100] axes, respectively. The reconstructed Ps D (P(oio)) below |p[ooi]| and 
b[ioo]| ~ 0.4 a.u. becomes positive at 12 K as seen in Fig. 3 (b). Highly positive values in 
the reconstructed p;? D (p(oio)) at 12 K appear in the region of |p(oio)| ~ 1-6 a.u. with peaks 
at |P[ooi]| = b[ioo]| ~ 11 a.u., similar to the situation at 300 K. 



IV. DISCUSSION 



In a 3d transition metal compound with spontaneous magnetization, the spin-magnetic 
moment is the difference in occupancy between the majority and minority spin bands, which 
mostly have 3d character. Electrons in fully occupied majority and minority spin bands 
do not contribute to the spin-magnetic moment. However, at any point in momentum 
space spin components of these band electrons are not completely compensated because 
hybridization effects depend on the spin state of the band electrons. These uncompensated 
spin components contribute to p s (p) [26]. Accordingly, p s (p) is the spin-uncompensated 
component of all occupied band electrons and the area under J ma g(p 2 ) is the spin-magnetic 
moment of the compound. 

Fe3C>4 is ferrimagnetic below 858 K where the magnetic moments of the Fe A and B sites 
in the inverse spinel structure are antiparallel, resulting in a saturation magnetic moment 
of about 4 p^/formula unit, which does not change drastically at T v 21]. The magnetic 
electrons are mostly 3d electrons from the Fe sites. Consequently, the extracted difference 
in the reconstructed Ps' D (p(oio)) results from a change in the band electrons with mainly 3c? 
character. 

Since the symmetry of Fes04 above T v is cubic, the negative values below |p[ooi] I and |p[ioo] I 
~ 0.4 a.u. in the reconstructed pg D (ptoio)) at 300 K stem from electrons in the minority spin 
bands below |p| ~ 0.7 a.u., a region in the first Brillouin zone of the inverse spinel structure. 
Accordingly, these electrons contribute to the electrical conductivity of Fe 3 4 because they 
have itinerant character. The reconstructed Ps D (P(oio)) below |p[ooi]| and |p[ioo]| ~ 0.4 a.u. 
is positive at 12 K and the positive values in the reconstructed Ps D (p(oio)) i n the region of 
| P(oio) I ~ 1-6 a.u at 12 K are smaller than those at 300 K. The population of electrons with 
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itinerant character decreases discontinuously at T v reflecting the remarkable decrease in the 
conductivity |3j. Therefore, the observed difference in the reconstructed Pg D (P(oio)) provides 
clear experimental evidence that electrons with |p| < 0.7 a.u. in the minority spin bands are 
localized below T v and that the electronic structure of Fe3C>4 with the inverse spinel structure 
is half-metallic in character. The MI transition is characterized by the localization of the 
electrons in the minority spin bands which contribute to the spin magnetic moment of Fes04. 
This MI transition with the structural change does not strongly affect the electrons with 
high momenta, which have a localized character even at 300 K. These results are consistent 
with a small change in the saturation magnetic moments at T v [2l|]. The behavior of the 
reconstructed Pg D (p(oio)) at 12 K cannot be understood using the simple ionic model with 
the inverse or normal spinel structure. 

We have evaluated the number of electrons, n e , whose wave functions localize from the 
positive region in the difference of reconstructed Ps D (P(oio)) between 12 K and 300 K. The 
value of n e is found to be 0.590(2) per unit cell (containing two Fe B sites) and should be 
compared to the charge ordering on the iron B sites of about 0.4 electrons per unit cell 
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16]. The higher value of n e indicates that the iron-oxygen hybridization plays 



also a role in the MI transition. 

We have carried out KKR band structure and momentum density calculations on Fe3C>4 in 



the inverse spinel structure within the framework o: 



the local density approximation (LDA) 



3. 



A reasonable level of agreement is 



in order to gain insight into the nature of p^ D (p) 
seen in Fig. 1 between the computed J ma g(Pz) and the measured spectra at 300 K along the 
three cubic axes, especially above p z ~ 2 a.u. Figure 3(c) shows the theoretical p1 D (j>) after 
it has been broadened to account for experimental resolution. The features of the calculated 
Ps D (P(oio)) look similar to those of the reconstructed Ps D (p(oio)) at 300 K. However, in 
contrast with the experimental spectrum, the positive peaks at |p[ooi]| — |P[ioo]| ~ 1-2 a.u. in 
Ps D (p) are sharper in the calculations and the value of the calculated Ps°(p) at |p(oio)| ~ 
is almost zero. These discrepancies cannot be understood within the ionic picture, but could 
arise from exaggerated hybridization in the LDA and/or correlation effects missing in our 
LDA-based calculations. As shown in Fig. 4, our computed spectrum before it is convoluted 
with the experimental resolution does display a negative value of the spin momentum density 
in the first Brillouin zone due to the presence of minority t<i 9 itinerant electrons appearing 
above T v , but this negative excursion is smeared out by the effect of resolution broadening. 
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The observed deeper crater-like feature in the reconstructed p;P(p(oio)) could be explained 
by the partial occupation of d 3z 2_ r 2 orbitals reflecting correlation effects beyond the LDA, 
since the majority spin d 3z 2_ r 2 orbitals give a positive contribution to Ps D (P(oio)) at |P(oio)| 
~ 0. Correlation effects are known to produce partial occupation of the natural orbitals and 
smearing of the electron momentum density 28]. 



V. CONCLUSION 



We have carried out magnetic Compton scattering experiments to investigate the spin- 
uncompensated electron momentum density p s (p) of Fe3C>4 at 12 K and 300 K. Mechanical 
de-twinning was used to overcome severe twinning in the single crystal sample below T Y . The 
reconstructed Ps D (P(oio)) below |p[ooi]| and |p[ioo]| ~ 0.4 a.u. changes from being negative 
at 300 K to positive at 12 K. This result reveals that electrons in the minority spin bands 
with |p| < 0.7 a.u. are localized below T v , while electrons with high momenta, which have 
a localized character at 300 K, are not affected strongly by the MI transition. Interestingly, 
the observed behavior of the reconstructed Pg D (p(oio)) at 12 K cannot be understood within 
the simple ionic model or the standard LDA-based band theory picture, and indicates the 
presence of correlation effects beyond the LDA leading to particularly strong smearing of 
the spin-uncompensated momentum density distribution. 
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FIG. 1: JmagiPz) of Fe304 along the three cubic principal axes as a function of p z . Open and filled 
circles with error bars represent the measured J ms ,g(p z )s at 300 and 12 K, respectively. Solid lines 
give the calculated J m ag(p?)s in the cubic spinel structure which have been broadened to reflect 



experimental resolution 



22J. 
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FIG. 2: Measured J ms , g (p z ) of Fe3C>4 between the cubic [001] and [101] axes in intervals of 9°: (a) 
at 300 K and (b) at 12 K with mechanical detwinning. 
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FIG. 3: (color online) (a) and (b) Reconstructed Ps (P(oio)) s of Fe304 from measured J ma , g (p z )s at 
300 K and at 12 K, respectively, (c) Calculated /fs Z) (P(oio)) of Fes04 in the inverse spinel structure 
with cubic symmetry which has been broadened to reflect experimental resolution. 
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FIG. 4: Reconstructed and calculated Ps D (P(oio)) s of Fe304 along the cubic [101] axis at 300 K and 
in the cubic spinel structure with and without convolution of experimental resolution, respectively. 
Filled circles represent the reconstructed Ps D (P(oio)) a * 300 K. Solid and broken lines give the 
calculated ^^(p^io)) 8 i n the cubic spinel structure without and with convolution of experimental 
resolution, respectively. 
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